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Abstract 

Vector and tensor analysing powers of the dp — ► {pp)n (charge-exchange) and 
dp — ► (pn)p (non-charge-exchange) breakup reactions have been measured with 
the ANKE spectrometer at the COSY ring at a deuteron beam energy of 1170 MeV 
for small momentum transfers to the low excitation energy (pp) or (pn) systems. A 
quantitative understanding of the values of A xx and A yy for the charge-exchange 
reaction is provided by impulse approximation calculations. The data suggest that 
spin-flip isospin-flip transitions, which dominate the charge-exchange breakup of 
the deuteron, are also important in the non-charge-exchange reaction. 
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It was suggested several years ago that the charge-exchange breakup of medium 
energy deuterons should show significant polarisation effects [1,2]. The tensor 

— * 

analysing powers of the dp — > (pp)n reaction were predicted to be especially 
large for excitation energies E pp of the final proton-proton pair below about 
3 MeV, provided that the momentum transfer q from the incident deuteron 
to the di-proton also remained small. Under these conditions, the final-state 
interaction in the ^So state of the pp system is very strong. The impulse ap- 
proximation predictions were tested successfully for very small E pp at beam 
energies of T d = 1.6 GeV and 2.0 GeV using both hydrogen and deuterium 
targets [3,4,5]. A refined numerical evaluation [6] proved that the theoretical 
description was equally valid over a larger E pp range at 200 and 350 MeV [7]. 
The cross section and tensor analysing powers are so large that the reaction 
could form the basis of a deuteron polarimeter with a high figure of merit [8], 
which was used to measure the polarisation of the recoil deuteron in large 
momentum transfer elastic electron-deuteron scattering at JLab [9]. 

One major feature of the reaction is that the differential cross section and 
the deuteron Cartesian analysing powers A xx and A yy are directly related 
to the magnitudes of the spin-spin neutron-proton charge-exchange ampli- 
tudes [2] . These govern the spin-transfer observables in the np — ► pn reaction 
at small momentum transfers from neutron to proton. This raises the possi- 
bility of using the deuteron charge-exchange reaction to obtain information 
about pn observables above say 800 MeV where, despite much work carried 
out at Saclay [10], accurate polarisation data are more sparse [11]. To lend cre- 
dence to such an approach, it is necessary to make detailed dp — > {pp)n mea- 
surements in a domain where one can have confidence in the nucleon-nucleon 
amplitudes used in the analysis. We therefore report here on the study of A y , 
A xx and A yy at T d = 1170 MeV (p d = 2400MeV/c) up to q = 130MeV/c. 
Quantitative agreement with predictions based upon an up-to-date phase 
shift analysis [11,12] is obtained, and this will allow the programme of using 
the charge-exchange reaction for pn studies to go ahead at COSY [13]. 

Though the situation regarding the dp — > (pp)n reaction looks fairly clear, that 
of the non-charge-exchange dp — > (pn)p breakup is far more complex. Even 
in the S-wave one has to consider the production of both and 3 S 1 (pn) 
states and in general the number of relevant low energy pn final states is twice 
that of pp, as is the number of elastic AW amplitudes that provide the driving 
force in the process. In a first exploration, we also took data in parallel on the 
dp — > (pn)p reaction in the range 70 < q < 200MeV/c, where the excitation 
energy E pn in the final pn system was similarly constrained to be at most a 
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few MeV. The results show some similarities to those of dp — > (pp)n under 
similar kinematic conditions and reinforce the belief that, at low momentum 
transfers, this reaction is also dominated by spin-flip isospin-flip transitions. 

The experiment was carried out at the COSY COoler SYnchrotron of the 
Forschungszentrum Jiilich using the ANKE magnetic spectrometer, located at 
an internal target position of the storage ring. Although ANKE has several de- 
tection possibilities [14], only those of the Forward Detector (FD) system were 
used here to detect the two fast protons from the dp — * (pp)n reaction [15]. The 
FD consists of multiwire chambers for track reconstruction and three layers of 
a scintillation hodoscope that permit time-of-flight and energy-loss determi- 
nations [16]. The tracking system gives a momentum resolution of better than 
1%. While fast protons from the dp — > (pn)p reaction were also measured 
in the FD, the slow recoil protons in the energy range 2.5 < T p < 32 MeV 
were detected in a Silicon Tracking Telescope (STT) placed inside the target 
chamber [15,17,18]. This provided a polar and azimuthal angular coverage of 
75.6° < # S tt < 116.4° and -21.8° < </> ST t < 19.4° respectively. 
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Fig. 1. Missing masses of the p(d,2p)X reaction obtained by measuring a) both 
final fast protons in the FD, and b) a fast proton in the FD and a slow one in the 
STT. Though there is slightly more background in case a), and the Gaussian peak 
is much wider (a = 13MeV/c 2 compared to 4MeV/c 2 ), there is no difficulty in 
extracting the p(d, 2p)n signal in either case. Events were retained in ±2.5 a regions 
around the peaks. 

In the deuteron charge-exchange reaction, two fast protons are emitted in a 
narrow forward cone with momenta around half that of the deuteron beam. 
As described in Ref. [15], such coincident pairs could be clearly identified us- 
ing information from the FD system, the coverage in the laboratory polar 
angle being between 0° and 6°. Because of this limited acceptance, there are 
strong kinematic correlations between the magnitudes of the two momenta for 
two- and three-body final states. For charge-exchange candidates selected in 
this way, the times of flight for the two particles from the target to the ho- 
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doscope were calculated, assuming the particles to be protons. The difference 
in these times of flight could be compared with the measured time difference 
for those events where the particles hit different counters in the hodoscope. 
This selection, which rejects about 20% of the events, essentially eliminates 
background, for example from dp pairs associated with n° production. The 
resulting missing-mass distribution for identified ppX events shows a clean 
neutron peak in Fig. la at M x = 940.4 ± 0.2MeV/c 2 (a = 13MeV/c 2 ), sit- 
ting on top of a slowly varying 2% background. 

The non-charge-exchange breakup dp — > {pn)p reaction was isolated by first 
identifying slow protons emerging from the target in the STT and then looking 
at the momenta of charged particles detected in coincidence in the FD. There 
is a gap of at least 800MeV/c between high momentum, elastically scattered, 
deuterons and protons from the dp — > (pn)p reaction, which have about half 
the beam momentum. The angular and energy resolution in the STT leads to 
a good excitation energy determination of the fast pn pair and, as a result, 
the missing-mass resolution on the neutron peak in Fig. lb is better than in 
la, with essentially no background. 

The COSY polarised ion source was set up to provide a sequence of an unpo- 
larised state, followed by seven combinations of deuteron vector (P z ) and ten- 
sor (P zz ) polarisations, where z is the quantisation axis in the source frame of 
reference. The determination of the actual polarisations of the beam, through 
the measurement of a variety of nuclear reactions, is detailed in Ref. [15]. It 
is shown there that the values of P z were on average about 74% of the ideal 
figures that could be obtained from the source, whereas the corresponding 
reduction factor for P zz was typically around 59%. 

Having identified the dp — > {pp)n events, these were binned in intervals of 
di-proton excitation energy E pp and three-momentum transfer q = and 
corrected for luminosity with the help of the beam current information in order 
to evaluate the analysing powers. In the right-handed coordinate system of 
the reaction frame, the beam defines the z-direction while the stable spin axis 
of the beam points along the y-direction, which is perpendicular to the COSY 
orbit. The numbers N(q, 0) of di-protons produced at momentum transfer q 
and azimuthal angle with respect to the x-direction are given in terms of 
the beam polarisations by 

N(q, 0) = No(q) [l + \P z A y {q) cos + \P zz {A yy {q) cos 2 + A xx (q) sin 2 0}] , 

(1) 

where N (q) are the numbers for an unpolarised beam, and A y (A yy , A xx ) are 
vector (tensor) analysing powers of the dp — > (pp)n reaction [19]. 

Measurements were made using the eight available spin states of the source. 
For each of the seven intervals in q (about 20MeV/c width), the yield was 
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Fig. 2. (^-distribution for the combination of counts (N 3 — N 2 ) (lower histogram) and 
(N 3 + N 2 ) /2 (upper) for spin states 2 and 3 for E pp < 1 MeV and 40 < q < 60 MeV/c 
of the di-protons from the dp — > (pp)n reaction at 1170 MeV. The crosses denote the 
calculated values of (AT 3 — _/V 2 ) from Eq. (1) using the known values of P zz in states 
2 and 3, the extracted tensor analysing powers, and the values of Nq = (iV~3 + N 2 )/2 
for the different bins. 

sorted in 120 bins of 3° width in 0, corresponding roughly to the resolution of 
the spectrometer. Using the known values of the beam polarisations P z and 
P zz [15], together with the average values of the trigonometric functions within 
each 0-bin, it is possible to determine the three observables A y , A xx , and A yy 
in a fit of Eq. (1). For each of the seven g-intervals, the resulting information 
from the eight different spin states was combined using a weighted average. 

The effect of the beam tensor polarisation can be seen immediately just by 
comparing spin states 2 and 3, where the ideal polarisations are (P z , P zz ) = 
(1/3, 1) and (—1/3,-1), respectively. The ^-distributions of two combinations 
of the number of events for these states are shown in Fig. 2 for E pp < 1 MeV 
and 40 < q < 60MeV/c. When events from all four quadrants, at 0, 180° — 
0, 180° + 0, and 360° — 0, are considered together, the term depending on 
A y is largely cancelled. In the ideal case, the average of the counts (N 3 + 
N 2 )/2 corresponds to that of an unpolarised beam, No, with a 0-dependence 
that reflects the acceptance of the spectrometer, independent of the beam 
polarisation. The difference (N 3 — N 2 ) is strongly modulated by polarisation 
effects. The crosses in Fig. 2 show the calculated values of (N3 — N 2 ) using the 
known values of P zz for the spin states 2 and 3, the extracted tensor analysing 
powers, and the measured values of (N 3 + N 2 )/2. 

In impulse approximation the vector analysing power is predicted to vanish 
for small E pp [2] and our results are consistent with this. The averages over 
the whole range in q are (A y ) = -0.001 ± 0.005 and (A y ) = -0.004 ± 0.004 
for 0.1 < E pp < 1 MeV and 1 < E pp < 3 MeV, respectively. The 0.1 MeV lower 
limit arises from the requirement that the protons hit different counters of the 
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FD. In order to improve the precision for the comparison with theory, we then 
imposed A y = and extracted values for the two tensor analysing powers. In 
this way the error bars are reduced compared to the fit with unconstrained 
A y , in particular for the highest q bins. The results are shown in Fig. 3 as 
functions of q, separately for the two E pp bins. The resolution in E pp is about 
0.3 Me V at 3MeV and better at lower E pp [20]. Due to the limited ANKE 
angular coverage, the acceptance gets steadily poorer as q and E pp increase, 
so that values of A yy could only be determined for q < 130MeV/c, with a 
slightly lower limit in the case of A xx . As explained later, the data are further 
divided into two groups depending on the angle 6 qk between q and the final 
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pp relative momentum k. 




Fig. 3. Cartesian tensor analysing powers A yy (open symbols) and A xx (closed) of 
the dp — »■ (pp)n reaction for a) 0.1 < E pp < 1 MeV, and b) 1 < E pp < 3MeV. 
The circles correspond to events where \cos6 q k\ < 0.5 whereas the stars denote 
| cos 6*gfc | > 0.5. The solid and broken curves, which involve respectively the same 
angular selection, follow from the impulse approximation program of Ref. [6], for 
which the 585 MeV input amplitudes were taken from Ref. [11,12]. The error bars 
include the uncertainties in the beam polarisation of about 4% [15]. 

Turning to the dp — > (pn)p reaction, the positioning of a relatively small aper- 
ture telescope in the horizontal plane means that only A y and A yy could be 
measured here. Furthermore, since (cos 2 0) ~ 0.96, there is a small contami- 
nation in A yy from the A xx contribution. In the extreme case that this is large 
and opposite in sign, as found for the charge exchange reaction in Fig. 3, this 
introduces a correction of about 6% to the values of A yy . These are shown in 
Fig. 4, along with those for A y , for events with E pn < 5 MeV. Although the 
resolution in momentum transfer was very good, the limited statistics pro- 
duced by the single telescope meant that the data were grouped in g-bins 
of different size (20MeV/c or 30MeV/c). The 2.5 MeV threshold for protons 
in the telescope implies that the lowest g-bin starts at 70MeV/c and so the 
overlap with the charge-exchange data set is not extensive. It is found that 
A yy is negative but, in contrast to the charge-exchange case, its magnitude de- 
creases with increasing q. Furthermore, whereas the charge-exchange A y was 
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consistent with zero for small E pp , here, though small, it increases steadily in 
magnitude with increasing q. 
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Fig. 4. Cartesian analysing powers ^ (crosses) and A yy (circles) for the dp — > (pn)p 
reaction for £7 pn < 5MeV in 20MeV/c or 30MeV/c bins in momentum transfer. 
The curves correspond to the predictions of A yy from the charge-exchange im- 
pulse approximation program of Ref. [6], as used for Fig. 3; E pn < 1 MeV (solid), 
E pn < 5 MeV (chain), and E pn < 9 MeV (dashed). 

The impulse approximation description of the cross section and analysing 
powers for the dp — > (pp)n reaction was developed in refs. [1,2,6] and we 
follow these works closely. To show the basic sensitivity of our measurements, 
consider neutron-proton charge-exchange amplitudes in the cm. system: 

fn P = a+i^(a- n +a- p )-h+P(a- n -h)(a p -h)+S(a- n -jh)(a p -jh)+e(a (2) 

where cr n and cr p are the Pauli matrices for the neutron and proton. The basis 
vectors are defined in terms of the initial (p) and final (p') momenta and lie 
along n = p x p' , I = p' + p, and m = n x I. 




For small momentum transfers and low excitation energy E pp of the final pp 
pair, the dp — > (pp)n charge exchange reaction mainly excites the ^So state of 
the final pp system. The spin-flip from the pn triplet to pp singlet provides a 
spin-filter mechanism. In single-scattering approximation, the resulting am- 
plitude depends upon the spin-dependent parts of f np , i.e. (3, 7, 5 and e but 
not the spin-independent term a. If, purely for the purposes of presentation, 
the deuteron D-state is neglected, at low E pp we expect that A y = and 

|/?| 2 + | 7 | 2 + k| 2 -2H 2 H 2 + | g | 2 -2|/3| 2 -2| 7 | 2 

l/3| 2 + |7| 2 + H 2 + kl 2 ' VV W + \l\ 2 + W + \e\ 2 ' 1 } 



Since (3 — 5 and 7 = at q = 0, the value of A xx = A yy depends there only on 
the ratio of \(3\ to which is believed to change smoothly with energy [2]. 



7 



However, the 8 amplitude, which contains the one-pion-exchange pole, varies 
very rapidly with momentum transfer and almost vanishes when q « m w . 
Hence A xx should approach its kinematical limit of +1 in this region, and this 
is consistent with the trend of the data shown in Fig. 3. 

Spin-triplet final states generally produce opposite signs for the analysing 
powers to the singlet states and so, away from the very small E pp region, there 
is dilution of the A xx and A yy signals. The estimation of this effect depends 
sensitively upon the pp final state interactions. In the program of Ref. [6], 
strong interactions [21] were kept only for L < 2, with the higher waves being 
distorted solely by a point Coulomb force. The deuteron S- and .D-state wave 
functions were taken from the Paris potential [22] but it was verified that the 
use of more modern potentials for the pp and np systems does not lead to any 
noticeable changes. The predictions were made using amplitudes derived from 
the current SAID NN phase shift solution [11,12]. 

It is known that, when q and the pp relative momentum k are perpendicular, 
odd partial waves cannot be excited and the pp system must be in a spin- 
singlet state [2] . As a consequence, less triplet dilution of the analysing powers 
is expected for small cos9 qk . To show this, we have divided the data shown in 
Fig. 3 into the two regions where | cos 9 q k\ §s 0.5 and imposed the same cuts 
on the theoretical description. All the features of both A xx and A yy are then 
reproduced, including the variation with q, E pp , and cos9 q k- It seems therefore 
that the model is as valid here as at lower energies [7] and that, as predicted in 
Ref. [2], multiple scatterings do not distort the analysing powers significantly. 

Returning to the dp — ► (pn)p reaction, the impulse approximation has not been 
evaluated in such detail for these data so that only qualitative statements can 
be made. If the deuteron D-state is neglected, the transition from the 3 Si 
deuteron bound state to the 3 <Si pn continuum is forbidden at q = due to the 
orthogonality of the two wave functions. The only allowed transition at q = is 
3 <Si — > ^o, which has a (AS, AI, AI Z ) = (1, 1, 0) character. It is the isobaric 
analogue of the deuteron charge-exchange amplitude and should exhibit the 
same analysing powers. Furthermore, the (0, 0, 0) transitions, driven by the 
large isoscalar spin-non-flip NN amplitudes, vanish like q 4 at small q. This 
is because they excite final 3 D\ or higher S'-waves that are orthogonal to the 
deuteron wave function. The behaviour can be illustrated particularly clearly 
if one sums over all excitation energies E pn to obtain the AS = closure sum 
rule [23,24]: 

+ \a pp - a pn \ 2 [1 - S(q)]} , 

\ v " / pd^(pn)p 

(4) 

where the a represent the spin-independent NN amplitudes. Expanding the 
S-wave deuteron elastic form factor S(q) in powers of q, we see that the 
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AI = 1 combination (a pp — a pn ) contributes to order q 2 whereas cancellations 
for AI = means that (a pp + a pn ) first contributes at order g 4 . 

To order q 2 the only source of dilution of the A yy signal due to spin-triplet 
final states arises from 3 S\ — > 3 -Po,i,2 transitions, which also involve an isospin 
flip. Given that the final transition to this order, 3 S\ — > l P\, corresponds to 
spin-singlet final states, we would expect that, for small momentum transfers, 
Ayy should behave much as for charge-exchange. 

The simple picture must be modified to include the deuteron D-state, though 
the basic suppression of the scalar-isoscalar transitions at small q remains. 
We therefore show in Fig. 4 the predictions of the charge-exchange impulse 
approximation for three upper limits on the value of E pn . Though at the lowest 
q the Ay y predictions are qualitatively right, the analysing power decreases at 
higher q while the predictions, even with a high 9 MeV cut-off, merely level 
off. The charge-exchange program prediction for A y reaches barely 10% of the 
dp — > (pn)p results shown in the figure, though it must be stressed that vector 
analysing powers are generally very sensitive to interference effects with small 
amplitudes. 

When the cos 6 q k cut is applied, very little change in the values of A yy is found 
in Fig. 4. This very different behaviour to that of charge exchange is due to 
the major dilution effect here being engendered by the I = 0, L = triplet 
final state, which is clearly unaffected by the suppression of odd partial waves 
by the angular selection. On the other hand, there is a modest increase in A y 
for the larger | cos 6^1, but any understanding of the significance of this will 
have to await a more complete theoretical investigation. 

In summary, we have measured the Cartesian analysing powers A y , A xx and 

— * 

Ayy for the dp — > (pp)n reaction at a beam energy of 1170 MeV, under kine- 
matic conditions where both the momentum transfer and the pp excitation 
energy are small, and found good agreement with the predictions of the im- 
pulse approximation. By detecting other calibration reactions simultaneously 
in ANKE, systematic effects arising from possible changes in luminosity with 
polarisation state are avoided [15]. The implications of acceptance cuts must 
be studied in detail before absolute cross sections can be obtained. Such infor- 
mation will be vital for isolating the AI = contribution by comparing the 
dp — > {pp)n and dp — > (pn)p cross sections. 

The use of a single silicon telescope placed in the horizontal plane meant that 
only A y and A yy could be determined for the dp — > (pn)p reaction though in 
the future, with a system of four larger telescopes [13], this limitation will be 
relaxed. The results are far more complex than for charge exchange and a more 
complete theoretical calculation is required to understand them quantitatively. 
It is important in such a model that the spin-triplet final states be treated 
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consistently for the bound and unbound systems, which means that the tensor 
force has to be included also when evaluating the pn continuum wave functions. 

Having shown that the analysing powers in the charge-exchange breakup re- 
action can be well described in impulse approximation, these measurements 
will be extended to higher energies in order to provide neutron-proton spin- 
dependent information in more barren regions, and may be subsequently used 
as a polarimetry standard for the COSY energy region. Although the maxi- 
mum deuteron energy available at COSY is only 1.15 GeV per nucleon, the 
energy range can be extended to almost 3 GeV by inverting the kinematics 
and using protons incident on a polarised deuterium target [13]. Furthermore, 
when both the proton and deuteron are polarised, information will also be 
gained on the relative phases of the pn charge-exchange amplitudes. 

We are grateful to R. Gebel, B. Lorentz, H. Rohdjefi, and D. Prasuhn for the 
reliable operation of COSY and the deuteron polarimeters. Neutron-proton 
scattering amplitudes corresponding to the current SAID phase shifts were 
kindly provided by I.I. Strakovsky and we are also indebted to R.A. Arndt for 
discussions regarding the SAID observables. 
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